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, The moments of the heavy quark-parton distribution functions in a heavy pseu- 

' doscalar meson, obtained in QCD sum rules, are expanded in the inverse heavy 
quark. Comparison with the finite mass results reveals that while the heavy mass 

r~| . expansion works reasonably well for the b quark, one has to take into account terms 

Qh| of higher than (l/rric)^ order for the c quark. 

1^ , In this talk I will very shortly report the result, obtained recently in my paper [1]. The 

problem is that there exist a large number of different theoretical models and approaches 

^ ' usually used to describe the heavy quark fragmentation functions (FF)(see, for example 

^ ' [2 — 6]), and expansions, based on the heavy-quark mass limit or on HQET are widely 

^ ■ used in many of them (e.g. [6j). So the question arise, does the expansion work good, 

T^ij- ! especially for charm case. In [1] it was offered to check the applicability (and accuracy) 

^ I of this expansion using the moments of the heavy-quark structure obtained long ago from 

On ' QCD sum rules in [7]. In this approach the moments of the heavy-quark parton distri- 

! 1 

, bution functions = J dxx^~^c{x, Q"^) were directly calculated in QCD sum rules and 

O ■ 



FF were estimated, following the familiar relation [S] ( see also [2] and [7] for details). 
The advantage of this method is the possibility to estimate the heavy parton distribution 
functions (hence, also fragmentation functions) in full QCD, in terms of universal param- 
eters, such as quark masses and condensate densities. One can use this result to fix the 
parameters of different models, as was discussed in [7]. 

There is no time to discuss the method itself, one can look it in the paper [7| and 
also [Ij. So I very briefly remind the main points of the method. One starts from the 
four-point correlator: 

U^, = tJ e*^^^'+*''^™^d^xdV'^(0 I T {j5(x)j;;^(y)jr^(0)jl(z)} I 0), (1) 

where = cF^c, and heavy-light pseudoscalar currents js = cT^u interpolate D meson. 
The above correlator is considered in the deep spacelike region at t = (pi — ^2)^ = 0, 
p^, qf < 0. We take into account the contributions of the unit operator (bare loop), quark 
condensate a = — (27r)'^(0 | ipip \ 0) and quark— gluon condensate (0 | ipG^^a'^^ip I 0). 
The gluon and four-quark condensate contributions are neglected, because the estimates 
show, that they are very small. We express this amplitude in terms of double dispersion 
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relation over and pi, then, us usual, saturate the amplitude by hadronic states and 
equate it to the result of OPE and perform two independent Borel transformations on pi 
and pI- (We will equate the Borel masses Mf^ = = only in the final sum rule). 
Following [7j, we choose the value of the variable s = (pi + gi)^ in the unphysical region 
as s = m|) — Q^, and transform the integrals into the moments of the parton distribution. 
Finally the result, obtained in, is [7j: 




1 + Z- m|5/g2)n+%2 + 4T^2^n+l 



Mn = Qd^D { dUi dz .^ , ° , . + , + ^6 (2) 
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Here 77 = (1 — (m^ — m\))/Q'^). The factor L accounts for the anomalous dimension of 
the quark condensate, normalized at point /x. We adopt the values of condensates P [TU]. 
normalizing them at /i = 1 GeV. In this case the factor L is close to the unity. 

In [7], it was shown that the first few moments of the sum rules are well behaved 
and the following numerical results were obtained at Ql = 20 GeV^ (Qg = 10 GeV^): 
M2 = 0.85(0.9), M3 = 0.75(0.83), M4 = 0.67(0.78) These moments can be used as an 
input in the evolution equations, at some initial point = {3 — 5)m|,, hence it is 
possible to predict the structure functions also at large Q^. 

Now we can answer to our question if it is possible to use the heavy— quark limit for 
the heavy quark parton distribution (and, correspondingly, fragmentation function). For 
that, we expand the sum rules for the moments in the inverse mass of the heavy quark. 

Substituting in Eq. (4) instead of the gj^imj^ the sum rule obtained from the two-point 
correlator [TT] we can rewrite the sum rules for Mn as 



(5) 



where the terms Rd{n) in the numerator and Kd in the denominator originate from the 
dimension c? = 0, 4, 6 contributions to the OPE of the four-point and two-point correlator, 
respectively. Following [7| we adopt the relation M^/2 = between the Borel masses 
in the four-point and two-point sum rules. Note that the exponential factor containing 
m\) cancels in the ratio of the two sum rules. 
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To investigate the heavy— quark mass hmit of ([5]) we employ the standard scahng 
relations for the heavy hadron mass, continuum threshold, and Borel parameter: rriD = 
m + A, So = TTi^ + 2muj, = 2rm and in addition assume the scaling = 'jmj^, 
where in the heavy— quark limit the parameters A, c<j, r and 7 do not depend of the heavy 
quark mass. (One should note, that A here has nothing common with well-known QCD 
parameter Aqcd)- In what follows we will use notation 6 = A/m, y = A/t.. 

Expanding in powers of 1/m , we obtain the first three terms for Rq: 



Ro{n) - D2 + (n + 1)5 



2D2 - (2 + 7)^3 
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^ ' 7^ + 47 + 3 + 



2{n + 2) + 37 
72(n + 2) 
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7V {n + 2)-fJ 372 ' {n + l){n + 2)^ 
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Here, D„ is defined as Dn = Eniuj/r)/ y"^^ and En{z) = J dze~^z^ For the contribution 



of (i = 4 operators we obtain: 



RAn) 



rL l + 2(5(n + l)/7 + (5' 



(n + l)(n + 2) 



[2/7 - 3/(n + 2)] . (7) 
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Finally, the contribution oi d = Q operators transforms into (first three terms) 



RM 



a m, 
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2S\n + l) 



{n + 2)[2 + {y + 4)/7 + 2/(^7)] - 3?/ - 24 - (8/3) (2n + ^)/y 
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And, finally, for the denominatori^ = Kq + + in ([5]) we obtain after expansion 



K c^Do- 2D^5 
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16r2 / 12r3 IQt^ y ' 



(9) 



From these equations it is easy to calculate the actual value of moments in the first 
and second orders of expansion and compare it with the exact value to examine the 
accuracy of the expansion in the inverse heavy quark mass. For numerical analysis we 
choose A = 0.6 GeV, u = 1.4 GeV, and r = 0.6 GeV, as usual. The numerical results 
for heavy— quark mass expansion of the first is shown in Fig. 1 where it is plotted as 
a function of the dimensionless ratio k = m/rric, for two different values of the ratio 
7 = Q'^/M'l) = 2,4. One can easily see, that even at relatively large masses of order 
of the 6— quark mass {k ~ 3) the heavy mass limit and even the first— order expansion 
in inverse mass are not reliable. Including the second 0{l/m'^) term one improves the 
situation for the 6— quark, but the approximate result for the moments still remains very 
far from the exact answer in the case of c— quark. We conclude that the heavy quark limit 
is not a reliable approximation for the parton distributions and fragmentation functions 
of c-quark. The analysis for higher moments (second and third) totally confirm this 
conclusion. 
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Figure 1: The first moment as a function of the quark mass m in the units k = m/rric. 
Curves 1,2,3 correspond to the first-and second-order expansion in the inverse quark mass 
and to the exact answer, respectively. The value of 7 = Q'^/M'^ = 2(4) is chosen in the 
left (right) panel. 
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